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Resonant exchange qubits are a promising addition to the family of experimentally implemented 
encodings of single qubits using semiconductor quantum dots. We have shown previously that it 
ought to be straightforward to perform a CPHASE gate between two resonant exchange qubits with 
a single exchange pulse. This approach uses energy gaps to suppress leakage rather than conventional 
pulse sequences. In this paper we present analysis and simulations of our proposed two-qubit gate 
subject to charge and Overhauser field noise at levels observed in current experiments. Our main 
result is that we expect implementations of our two-qubit gate to achieve high fidelities, with errors 
at the percent level and gate times comparable to single-qubit operations. As such, exchange-coupled 
resonant exchange qubits remain an attractive approach for quantum computing. 


I. INTRODUCTION 

The seminal work of Loss and Divincenzo^ introduced 
the notion of using individual electrons trapped in gate- 
defined quantum dots to encode quantum information, an 
idea which has since burgeoned into a family of promising 
architectures for quantum computing [Mg. 

An early theoretical realisation was that a single qubit 
encoded in three electron spins could be universally con¬ 
trolled using exchange couplings alone m, which re¬ 
moves any requirement for individually addressable elec¬ 
tron spin resonance or magnetic field gradients. Cru¬ 
cially, since exchange couplings in semiconductor experi¬ 
ments are controllable using gate voltages, this allows all 
qubit operations to be performed electronically; an at¬ 
tractive feature in experimental implementations. Single¬ 
qubit operations for the so-called “exchange-only” qubit 
have been experimentally demonstrated nmn]. Pulse 
sequences are known for single-qubit gates that simul¬ 
taneously correct leakage errors and other sources of 
noise mi. Proposed two-qubit gates for the exchange- 
only qubit either involve capacitive coupling [isiin] or 
exchange-coupling [13 HMD]. Since exchange coupling 
is usually much larger than capacitive coupling, exchange 
gates are usually faster, but come at the cost of requir¬ 
ing complicated pulse sequences in order to echo away 
the unwanted spin-flip transitions that occur as a side- 
effect and cause leakage errors. The first such pulse se¬ 
quence m that effected a CNOT required 19 exchange 
pulses in 13 timesteps, and was found using a numerical 
search. Since then, robust numerical searches have found 
improved pulse sequences that are robust against more 
sources of decoherence and/or reduce the number of gate 
operations HSHiD]. 

A recent alternative to the “exchange-only” qubit is 
the “resonant-exchange qubit” mi ED, which encodes 
qubits in the interaction picture with respect to signi¬ 
ficant exchange coupling between the three dots. Uni¬ 
versal single qubit operations are effected using rf gate 
pulses to the electrodes controlling the exchange coup¬ 
lings. This qubit has been shown both theoretically m 


and experimentally mi to have several improved prop¬ 
erties, including first-order insensitivity to charge fluc¬ 
tuations and reduced leakage error due to nuclear field 
fluctuations III1E2]. In addition, Taylor et al. m have 
shown that you can perform two-qubit gates between 
these qubits using charge dipole interactions, while we 
have suggested an alternative of using simple exchange 
pulses between nearby qubits . 

In this earlier work [23] . we showed that a two-qubit 
CPHASE gate can be implemented using a single ex¬ 
change pulse between the constituent quantum dots of 
neighbouring qubits (shown schematically in figure |^. 
Rapid high-fidelity gate operation is in principle made 
possible by energetically suppressing the spin-flip trans¬ 
itions that lead to leakage. This method of effecting a 
two-qubit gate contrasted with the more conventional ap¬ 
proach of long and complicated pulse sequences, and is 
similar to our earlier proposal for singlet-triplet qubits 

E3]- 

It is the purpose of this work to extend our previous 
results by considering higher-order analysis, adiabatic 
pulse profiles, the effect of noise, and the performance 
of our gate in physically motivated simulations that in¬ 
clude noise. In section]^ we briefly review single-qubit 
resonant-exchange qubit operations, in section jllTI we re¬ 
view our two-qubit gates, in section [IV] we formally char¬ 
acterise our two-qubit gate in the butterfly geometry, in 
section [V| we briefly consider other geometries, and in 
section IVll we conclude. 

In this work we set h = 1, meaning that energies are 
interchangeable with angular frequencies. 

II. RESONANT EXCHANGE QUBITS 

A resonant exchange qubit is a triple-dot system op¬ 
erating deep in the (1,1,1) charge state (each dot al¬ 
most surely confines a single electron). A large magnetic 
field is applied along the z-axis, Zeeman splitting the 
23 = 8 spin states according to the z-projection of their 
total spins: AE = B± is the effective Zeeman 

splitting of an electron subject to the global magnetic 
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Figure 1 . (colour online) A Bloch sphere schematic depicting 
the effect of intra-qubit exchange couplings on qubit states. 
Under J12 and J23 couplings, states rotate about the same 
axes as those of the rotations indicated in green and red re¬ 
spectively. If both J12 and J23 are equal to J^, states rotate 
about the z axis as indicated by the emboldened blue rota¬ 
tion. Pulses Involving different J12 and J23 couplings allow 
arbitrary rotations, and hence provide universal single qubit 
control. As these intra-qubit couplings are controllable using 
gate voltages, this allows for universal electronic control of 
qubits. 


field. Note that we have absorbed the g-factor (up to 
its sign) into our definition of B±. The logical states 
of the qubit are | 0 ) = (Itti) + lift) ~ 2 /v^ and 

| 1 ) = (Itti) ~ litt)) /v^, which both have total spin 
S' = 1/2 with z-projection of ruz = 1/2. The remainder of 
the eight-dimensional Hilbert space describes non-logical 
states, the spanning eigenstates of which are completed 
for the energy eigenbasis in table These qubits are op¬ 
erated with large intra-qubit exchange couplings J 12 and 
J 23 , with oscillatory modulations around J 12 = J 23 = Jz 
providing single qubit control HU HU. This is depic¬ 
ted and described in more detail in figure Intra¬ 
qubit exchange couplings are in turn controlled by de¬ 
tuning the voltages defining the quantum dots (which we 
parametrise as e). This allows for complete electronic 
control of qubits, which potentially simplifies experi¬ 
mental implementation. In several GaAs singlet-triplet 
qubit experiments H El H Hg, an exponential ansatz 
J(e) = Joexp(e/eo) has been found to be a good phe¬ 
nomenological fit to experimental data over a wide range 
of interesting values of e, and so we adopt it in this work. 


III. EXCHANGE-COUPLED TWO-QUBIT GATE 

Here we provide a brief review of the two-qubit gate 
between resonant exchange qubits described in our prior 
letter [23]. Consider two resonant-exchange qubits {A 
and B) in a large transverse magnetic field, coupled in 
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State 

s 
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Bx/ 2 - 3 J ./2 

|i-) 

lut) - im) 

1/2 

-1/2 

B_l/2-Jz/2 

\Q-) 

| 4 Tt) + lltl) + ItlT) 

3/2 

-1/2 

Bx /2 

Q-3/2) 

lUl) 

3/2 

-3/2 

3 Bx /2 


Table I. Energy eigenstates of a single resonant exchange 
qubit, sorted by m^. S is the total angular momentum 
quantum number of the three electron spins and rriz is the 
z-component of the total angular momentum. Jz = J12 = J23 
is the energy level splitting due to exchange coupling, and 
is the splitting of an electron subject to the large transverse 
global field. Each eigenstate with a negative subscript has 
their constituent spins flipped relative to the corresponding 
unsubscripted state. 
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Figure 2 . (colour online) The physical arrangements of 
quantum dots (or geometries) considered in this work. Dashed 
green lines indicate large intra-qubit exchange couplings, 
with solid blue lines indicating the weaker exchange coupling 
between qubits that can be used to effect two-qubit gates. 
We will focus mainly on the (a) butterfly and (b) linear and 
geometries, with the (c) rectangular geometry being another 
alternative. 


several different ways as depicted in figure . We label 
the intra-qubit couplings of the qubits and , and 
assume that the couplings within each qubit are equal 
(i.e. The logical states of the two 

qubit system, | 0 , 0 ), | 0 , 1 ), | 1 , 0 ) and | 1 , 1 ), are all in the 
niz = 1 subspace, along with eleven other states. Since 
we will only be considering dynamics which conserve z- 
projection of spin rUz and are working in a large global 
magnetic field, this reduces the dimension of the Hilbert 
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Geometry SJzjJc JzzjJc. J^/Jc 


Linear 

1/36 

1/36 

- 1/24 

Butterfly 

- 1/18 

1/9 

0 

Rectangular 

0 

1/6 

-1/12 


Table II. Qubit coupling parameters arising from lowest order 

g srturbation theory in each of the three geometries of figure 
Calculations assume that j/ ~ J^, and that all non-zero 
exchange couplings Jij are equal to Jc. When \ — J^| 3 > 

Jc, the degeneracy of the logical |10) and |01) states is broken, 
and we find that Jx —>■ 0 for all geometries. All other entries 
in the table are unaffected. 

space that can interact with the logical subspace to a 
maximum of 15; though in some cases (such as for the 
butterfly geometry) symmetry constrains this subspace 
further. For an explicit representation of all 15 energy 
levels in the energy eigenbasis, please refer to the supple¬ 
mentary material. The most important observation to 
make about the non-logical = 1 subspace is that they 
all have energies different from those on the logical sub¬ 
space by at least min( , Sj_). This guarantees that 

leakage transitions will be unfavourable provided that all 
additional energy-level splittings remain less than ~ Jz, 
which leads to times of order at least ~ l/J^. Since 
single qubit operations already run slow compared to 
these timescales , this is not a restriction in practice. 

Using the fact that leakage processes are suppressed 
when Jc <C Jz, we performed first order degenerate 
perturbation theory around Jc = 0 and wrote down 
its effect on the logical subspace; the so-called effective 
Hamiltonian[^ on the logical subspace. The zeroth or¬ 
der terms arising from perturbation theory describe the 
uncoupled resonant exchange qubits. Writing the Pauli-Z 
logical operators on qubit A and B as af and af respect¬ 
ively, the zeroth order effective Hamiltonian is: 

iJo = -{Jt + Jz) - \jt<^z - \jz^z- (1) 

The first order terms describe the leading order effect 
of the inter-qubit coupling Jc- The effect on the logical 
subspace is described by: 

He = SJc + ^SJziaf + erf) 

+Jzzcrfaf + j±{a^a^ + a^a^), ( 2 ) 

with Jzz, SJz and J± all being geometry dependent, as 
specified in table |TTj This perturbative analysis is re¬ 
peated in greater detail in the supplementary material. 

The structure of these effective Hamiltonians admit a 
straightforward two-qubit CPHASE gate using a single 
DC exchange pulse. Simple AC coupling pulses may also 
be interesting |^, but we leave this to future work. Con¬ 
sider first two qubits coupled according to the butterfly 
geometry of figure |2b| Since J_l = 0 in this geometry, the 
only two-qubit component of the gate’s operation at first 
order is Jzz, which will implement a CPHASE gate after 


a time r such that fj Jzzdt = 7 r /4 (modulo single qubit 
unitaries). The other two geometries, linear and rectan¬ 
gular in figures [ 2 a| andrespectively, have non-zero J^\ 
and consequently will not perform a CPHASE gate unless 
the additional contribution can be suppressed. This can 
be achieved by detuning the intra-qubit exchange coup¬ 
ling energies such that | jf — jf | ^ J^; or by adding 
a simple logical Z (uf or uf) echo pulse at t = r/2 to 
one of the qubits associated with each exchange coupling 
(which anti-commutes with erf erf and erfaf, and thus 
cancels out the effect of J_l). We opt not to consider more 
sophisticated pulse sequences that echo out higher order 
contributions to Jx in order to maintain the simplicity 
of our gate. 

IV. GATE CHARACTERISATION 

We now begin a more complete characterisation of the 
performance of our two-qubit gate. The gate has two in¬ 
trinsic sources of error (which would be present even in 
an ideal implementation): timing inaccuracies and leak¬ 
age; and we will consider the two sources of extrinsic 
noise anticipated to be most pertinent in experimental 
implementation: charge and Overhauser noise. 

In the following two subsections, we will show that 
intrinsic noise can be effectively mitigated by correctly 
tuning gate times and by adiabatic pulse sequences, res¬ 
ulting in high fidelity gate operations. We then move 
on to consider how robust our gate is to the anticip¬ 
ated sources of experimental noise. Simulations will 
be provided in each section to demonstrate the anticip¬ 
ated performance of our gate. These simulations involve 
monte-carlo averaging (over pseudo-static parameters) of 
solutions to a lindblad master equation (encoding high 
frequency noise). The performance measure used is “en¬ 
tanglement fidelity”, as described in a former work m- 
“Entanglement fidelity” is related to the more commonly 
used “average fidelity” of random benchmarking by: 

- dFf. + 1 

where F is the average fidelity, F^ is the entanglement 
fidelity, and d is the dimension of the quantum system 
[271 [2E] (J = 4 for our two-qubit system). Entanglement 
fidelity is used in this work because it can be directly 
computed using a fixed input state, which simplifies sim¬ 
ulations. In this section, the simulations are usually done 
for the butterfly geometry which has the greatest sym¬ 
metry and performance. In the next section ([V|), we ex¬ 
tend our analysis to the linear geometry, which should be 
easier to fabricate for experiment. 

A. Timing Inaccuracies 

Due to the complexity of the dynamics of the six 
quantum dot system, there is no closed analytic form for 
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the ideal gate time. As a result, one needs to be careful 
how the gate time is estimated; over- or under-estimating 
the ideal gate time will result in a corresponding over or 
under accrual of two qubit phase, and thus reduced gate 
fidelities. 

The ideal gate time is the time r for which a noise¬ 
less exchange pulse should be turned on between the two 
triple-quantum-dot systems in order to perform a two- 
qubit CPH ASE gate on the encoded qubits. Recall from 


section III that r is implicitly defined by Jzzdt = 7r/4, 
with Jzz being geometry (and potentially time) depend¬ 
ent, as shown in table 

Using the butterfly configuration as an example, first 
order perturbation theory predicts that Jzz{Jc) = Tc/9 
and thus implies that r = Ott (Jc)^ /4; where {Jc)t ^^e 
time average of during the pulse. Note that the lin¬ 
earity of the first order approximation for Jzz(Jc) allows 
one to compute gate time r in a manner agnostic to the 
details of the pulse shape, requiring knowledge only of 
the average value of Jc during the pulse. This property 
is lost beyond first order, as Jzz{Jc) has corrections at 
higher order that become significant in all geometries for 
physically relevant values of Jc and Jz, meaning that r 
must be calibrated anew for each pulse shape. 

There is no closed analytic form for r, and since r 
would in any case have to be calibrated in-situ in any ex¬ 
perimental implementation using one of several optimisa¬ 
tion techniques |29f[52] , we refer the reader to the supple¬ 
mentary material for a description of how we numerically 
optimise r in our simulations. Henceforth, we assume r 
has been estimated perfectly, and note that in many of 
our simulations r differs significantly from its first order 
estimates. 


B. Leakage 

Leakage is a measure of how much a state initially 
with support only on the logical subspace shifts sup¬ 
port onto the non-logical subspace during a logical op¬ 
eration, and results in reduced entanglement fidelities 
(to first non-trivial order, J- = 1 — C). For an arbit¬ 
rary state p, we quantify this using C = Tr(PpP), where 
P = 1 —0110 11 10 (0 is the projector off the logical 
subspace. Leakage occurs via energetically forbidden ex¬ 
citations that are suppressed by the energy gap between 
the logical states and a leakage state, or when logical 
states are subjected to pulses with frequencies corres¬ 
ponding to the energy gap. 

In the two-triple-quantum-dot system, logical states 
are isolated from leakage states by an energy gap AP 
proportional to Jz- When coupled using Jc, this en¬ 
ergy gap monotonically reduces. Consequently, as the 
ratio Jc! Jz increases, the likelihood of leakage also in¬ 
creases. The energy level spectrum, with an indication as 
to which non-logical states the logical states can couple, 
is shown in figure Note that the symmetries of the 
butterfly geometry cause Jc to couple logical states to 
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Figure 3 . (colour online) The energy level spectrum of the 
niz — 1 subspace of a two resonant exchange qubit system. 
Energy levels are each labelled, with starred levels indicat¬ 
ing logical eigenstates. The eigenstates corresponding to each 
label are explicitly written out in the supplementary mater¬ 
ial. Under butterfly coupling Jc, the logical states couple to 
like coloured/dashed states leading to a minimal energy gap 
of ~ 3 Jz/ 2 . Note that the like coloured/dashed levels form 
distinct subspaces. Under linear coupling, all coloured/thick 
energy levels are coupled (including logical states), leading to 
a minimal energy gap of ~ Jz/ 2 . 


disjoint subspaces and give rise to an effective energy gap 
of AE = 3Jzj‘2, which is three times larger than in the 
linear system where AE = Jzl‘2- due to Jc coupling all 
of the logical states into the same subspace. This leads 
to substantially improved performance in the butterfly 
configuration for any given JdJz- 

Implementing our two-qubit gate requires Jc to be act¬ 
ive only for a fixed duration r, meaning that rapid or 
broadband changes in Jc when it is turned on and off 
can lead to excitations from the logical subspace. Choos¬ 
ing pulse shapes with discontinuities only at high differ¬ 
ential orders can therefore further suppress leakage by 
several orders of magnitude, and hence improve fidelit¬ 
ies, as shown in figure For more intuition regarding 
adiabatic pulses and leakage refer to our earlier work on 
adiabatic pulses for singlet-triplet qubits |24j . 

In the absence of noise, the entanglement fidelity J- « 
1 — A is limited only by leakage, which is in turn de¬ 
termined by the choice of adiabatic profile and the ratio 
Jc! Jz- This provides an upper bound on the performance 
of implementations of our gate. In this work, we will use 
a narrow-band sinusoidal adiabatic pulse described by: 


= Jr ( 1 - cos ( 


which has its first discontinuity at second order when 
t = 0 and T. Note that Jc is chosen such that Jc is the av¬ 
erage value of the pulse. This choice allows for fast single¬ 
pulse gates with fidelities in excess of 0.9999 for physic¬ 
ally reasonable parameters, as shown in figure and 
hence this gate may prove to be useful for fault-tolerant 
computation using semiconductor quantum dots. The 
remainder of this section will be devoted to determining 
how robust this performance is to anticipated sources of 
experimental noise. 
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Figure 4. (colour online) (a) Leakage L and (b) infidelity \—T 
(bottom) at the end of a single two-qubit gate operation for 
several different adiabatic profiles in both the butterfly and 
linear geometries. The x-axis is shared between the plots, 
and is over Jc/Ai?: the ratio of inter-qubit coupling and the 
geometry-dependent minimum energy gap. The upper x-axis 
provides a conversion from Jc/AE to Jd Jz for the butterfly 
geometry (an expression in terms of the controllable paramet¬ 
ers of the model). Colours indicate the adiabatic pulse pro¬ 
file used, while solid (dashed) lines denote that the butterfly 
(linear) geometry is being considered. Crucially, these plots 
demonstrate that use of adiabatic pulses can improve suppres¬ 
sion of leakage by several orders of magnitude, provided that 
Jc is small compared to the energy gap; and that this leads 
to a corresponding increase in gate fidelities. 


C. Charge Noise 


One of the most significant sources of experimental 
noise affecting semiconductor quantum dot qubits is 
charge noise mm- Charge noise is the effect of fluc¬ 
tuations in electric potential on the gates defining the 
quantum dots, for which there are many causes includ¬ 
ing environmental rf radiation and Johnson noise |35| . 
Charge noise on the electrodes exhibits itself in our model 
as fluctuations in the electrode detunings e. 

In this work, we describe charge noise on e using a 
two-parameter phenomenological model that approxim¬ 
ates the noise in gate voltages by static and white noise 
perturbations around the desired value. The pseduo- 
static (DC) and white noise (HF) perturbations are re¬ 
spectively parameterised by the standard deviation of the 
pseudo-static charge offset and the spectral density 


Parameter Value 


Exchange couplings: 

J, = = Jf 1.65 peV (~ 0.4 GHz) 

Exponential ansatz: 

Jo 82.7peV (~ 20GHz) 

Ed 0.35 mV 

Gharge noise: 

(Te 15.8 pV 

D 0.244 pV^ ns 

Overhauser noise: 




2mT 


Table III. Model parameters. In the experiments of Medford 
et al. m, the intra-qubit couplings Jz had measured val¬ 
ues of roughly 0.8peV (0.2 GHz) to 4peV (IGHz), of which 
we’ve chosen a conservative value. The parameters for the ex¬ 
ponential ansatz J(e) = Joexp(e/eo) were chosen to roughly 
match the experimental results of Dial and collaborators |25| . 
The noise parameters and D were calibrated by respect¬ 
ively matching somewhat typical values of T 2 = 25 ns and 
T 2 — 20ps from experiment |lllI15| . The standard deviation 
of the Overhauser field cts was chosen to be consistent with 

0 a El EH- 


of charge fluctuations D; which can be respectively cal¬ 
ibrated to experimental T 2 and T2 characteristic times. 
This is the same model described in our earlier work on 
singlet-triplet qubits El], which seems to reasonably de¬ 
scribe the results of experiments ESI, even though the 
precise mechanisms that cause this behaviour are not 
perfectly understood |35| . In experiment, there are likely 
to be additional high-frequency Ti processes biased to¬ 
ward relaxation, which we have chosen not to include 
in this model as it would require adding parameters to 
our model that have not been sufficiently well empirically 
constrained. Moreover, it has been found in experiment 
that characteristic Ti times are in excess of ~ 40 ps for 
Jz < 1.5peV m, which is long compared to T2 times 
of T2 ~ 20 ps m- Even though Ti times are found to 
decrease with larger Jz E], and so at some point will be¬ 
come comparable to the dominant sources of noise, the 
performance of our gate does not appear in any case to be 
limited by high frequency noise in the parameter space of 
interest to us (see figure |5b| ) . We therefore do not expect 
that the fidelity of our gate operations will be strongly 
affected by Ti relaxations. 

During the operation of our gate, there are multiple 
exchange couplings active at once (five in the case of the 
butterfly and linear geometries). We assume that charge 
noise is independent on each coupling, and use as men¬ 
tioned earlier an exponential ansatz for each coupling: 
Jij{sij) = Jq e'xp(sij/eo), where i and j indicate the pair 
of dots being coupled. We also assume that the size and 
spacing of each pair of dots is the same, allowing us to use 
the same Jq and Eq for each coupling. Jq and Eq are calib¬ 
rated to match such that J{e) is consistent with the data 
from Dial et al. E3. In the small noise limit in which we 
are interested, first order analysis of the effect of perturb- 
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ations to the intra-qubit couplings J12 and J23 allow one 
to derive that Og = JjT^ and D = jT^'^ 

where J is the exchange coupling inferred from the ex¬ 
periment. In this work, we fit these parameters to ex¬ 
perimental single resonant-exchange qubit T 2 times of 
^ 20 ps and T 2 times of ^ 25 ns [TTl [15]. The resulting 
parameters are shown in table jlllj 

The simplicity of the noise model and exponential an- 
satz allows us to construct a qualitative model from per¬ 
turbation theory of the effects of this noise model on the 
entanglement fidelity at the end of a single two-qubit gate 
operation: 


+ y/2Jl/Jl) {al + DJjk) ,(3) 

where C is the ultimate leakage, Jc and are the inter- 
and intra-qubit exchange couplings respectively, eo is 
parameter of the exponential ansatz and fc is a geometry 
dependent constant term. 

This qualitative model provides several key insights. 
Firstly, we learn from the (l-I-factor that 
the fidelity dimunition due to noise becomes more signi¬ 
ficant for longer gate times {Jc/Jz small), and so one is 
encouraged to run the gate as quickly as possible. The 
\/2 j term arises from DC noise on intra-qubit coup¬ 
lings, and proves to be the dominant contribution to infi¬ 
delity. The (cr^ -I- I? Jc/fc) factor suggests that DC noise 
(parameterised by a^) leads to an appoximately uniform 
dimunition of fidelity for any given Jc/Jz, whereas high 
frequency noise (parameterised by D) becomes relatively 
more significant as Jc increases; and that their effect is 
additive. Putting this all together we predict that gate 
performance decreases when Jc/Jz is too large (where 
leakage errors dominate) or too small (where low fre¬ 
quency charge noise dominates). For the experimentally 
relevant parameters considered in this work, we will not 
enter a regime where Jc is sufficiently large that high 
frequency charge noise dominates. These predictions are 
corroborated in simulations of the butterfly geometry in¬ 
cluding high frequency and DC charge noise, as shown in 
figure]^ With the noise parameters chosen to correspond 
to current experiments, we note that the gate’s perform¬ 
ance in the presence of charge noise appears to be limited 
by intra-qubit DC noise. This pseudo-static noise can in 
principle be echoed out by a single echo pulse. We will 
consider the effect of such simple echo pulses later when 
discussing the linear geometry. As well as echo pulses, we 
expect that technical developments will reduce the level 
of low frequency charge noise in the future. 


D. Overhauser Noise 

Perhaps the most widely used substrate in semicon¬ 
ductor quantum dot experiments is the GaAs/AlGaAs 
heterostructure. The nuclei in this substrate have non¬ 
zero spin, which give rise to a net magnetisation which 


slowly (compared to gate times) varies due to nuclear 
spin flip-flop interactions |3h| . This net polarisation is 
called an Overhauser field. While there exist semicon¬ 
ductor substrates composed of nuclei which do not have 
a spin (such as silicon and graphene), GaAs/AlGaAs has 
remained a popular material due to the well-developed 
fabrication techniques associated with it. In this sec¬ 
tion, we quantify the effect of varying levels of Overhauser 
noise on the performance of our gate. 

The Overhauser field looks like an additional randomly 
oriented local magnetic field at each dot. We assume that 
only a small fraction of the ~ 10® nuclei in the vicinity of 
each dot also contribute significantly to the polarisation 
of an adjacent dot, and therefore make the approximation 
that these local fields are uncorrelated. In GaAs/AlGaAs 
structures, the Overhauser field has an RMS magnitude 
of about 1 to 3mT m- 

We model the Overhauser field as a pseudo-static off¬ 
set of the magnetic field along the z-axis sampled from 
a normal distribution with standard deviation as- The 
neglected in-plane components of the random field con¬ 
tribute at second order in perturbation theory, and their 
effect is suppressed when the system is subject to a large 
magnetic field along z (as already posited in earlier sec¬ 
tions), and so can be safely ignored in this analysis. For 
a more complete analysis, we refer the reader to Hung et 

al. |57] . 

The effect of this random Overhauser field is to couple 
the logical subspace to all of the fifteen niz = 1 states 
shown in figure leading to both leakage and logical 
errors. The logical errors directly caused by the Over¬ 
hauser field are limited to single qubit errors, as must be 
the case since the magnetic field perturbations are local 
to each dot. 


Simulations of entanglement fidelity for a range of dif¬ 
ferent Overhauser field magnitudes in an otherwise noise¬ 
less gate implementation are shown in figure Our 
monte-carlo simulations preferentially sample some mag¬ 
netic field gradients in order to increase the rate of con¬ 
vergence, as described in the supplementary material. 
From we see that the effect of the Overhauser field 
is approximately linear in gate time (breaking down only 
once leakage becomes significant). From 6b we see its 
effect is also approximately linear in the standard devi¬ 
ation of the field. 


It should be noted that various techniques exist 
for reducing the magnitude of the Overhauser field in 
GaAs/AlGaAs, for example by nuclear state preparation 
which has been shown to reduce the RMS by a factor of 
~ 70 [3S] ; in which case the effects of the Overhauser field 
can be largely neglected. Implementing these techniques 
can be quite complicated, however, and so we assume 
conservatively that the field will be unsuppressed . 
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Figure 5. (colour online) The contribution to fidelity (T^noiay — T^noiseiess) from (a) pseudo-static (DC) and (b) white (HF) 
charge noise on each of the five exchange couplings active during a single gate operation using the butterfly geometry. The 
colourmap used is divergent at a fidelity of 0.9, as shown above; with blue colours indicating performance in excess of 0.9, and 
grey colors indicating performance below 0.9. Since DC noise is essentially dependent only on the the ratio Jc/Jz, we plot in 
(a) the infidelity contribution as a function of Jc/Jz and The value of from table III is indicated by a dashed line at 


15.8 In (b) we plot the fidelity contribution from white charge noise for D = 0.244 fiV^ns as in table 
behaviour of both plots is predicted by equation 


III 


The qualitative 



Ut^eV) 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 

(Tb {mT) 


Figure 6. (colour online) The contribution to fidelity (J^noisy — J^noiseiess) from a pseudo-static Overhauser field with standard 
deviation as. In (a) the fidelity contribution is plotted as a function of Jc/Jz and Jz for as = 2mT, revealing that Overhauser 
noise contributes to infidelity in a manner somewhat proportional to gate time (contours of constant gate time are indicated 
by grey lines), except at large Jc/Jz corresponding to large leakage. In (b) the fidelity contribution is plotted as a function of 
Jc/Jz (with fixed Jz = 1.65/rel^) and as. The radial contours of (b) (in regions where leakage does not dominate) imply that 
infidelity is proportional to as/Jc which is roughly oc ras, thus corroborating that the infidelity grows roughly as gate time. 
The colourmap used is the same as in figure 


E. Cumulative Noise Model 

As mentioned earlier, charge and Overhauser noise are 
the two most significant sources of noise in semiconductor 
quantum dot experiments, and in this section we simulate 
the performance of our gate in the presence of both. 

In principle the effects of Overhauser and charge noise 
can be more serious than for either noise source alone. 
For example, we learn from perturbation theory that 


at second order and above, the Overhauser field cross¬ 
couples with charge noise in Jc, allowing the Overhauser 
field to effect genuine two-qubit errors. The larger J^,/ Jz 
and (Tb become, the greater this cross-coupling and hence 
the non-linearity in the contribution of charge noise and 
Overhauser noise to fidelity. Fortunately, this effect is 
small in the parameter regime in which we are interested, 
and so we refer the reader to the supplementary inform¬ 
ation for more details. This is visible in the simulations 




















































perimental evidence m suggests that the parameters 
we have chosen are not too far from optimal for cur¬ 
rent devices. Secondly, one could add a spin echo pulse 
at t/2, which will be discussed in more detail for the 
linear geometry in the next section. For the butterfly 
geometry, simulations involving an echo pulse boost fi¬ 
delities to ~ 0.98 at the cost of slowing down the gate. 
Thirdly, fidelities could be increased by suppressing the 
Overhauser held using DNP, or using a material without 
nuclear spin. We found that suppressing the Overhauser 
held entirely while also echoing the low frequency charge 
noise allowed for fidelities exceeding 0.99. 


Figure 7. (colour online) Entanglement fidelity for the but¬ 
terfly geometry subject to charge and Overhauser noise, as 
described in the main text, with noise levels calibrated to cor¬ 
respond to experiment (see table III I. Performance is found 
to be limited predominantly by leakage (from above), intra¬ 
qubit DC charge noise (from below), and Overhauser noise 
(from the left); leading to an optimal ratio of Jc/Jz ~ 0.6 
where fidelity reaches ~ 95% for large enough Jz. The col- 
ourmap used is the same as in flgure[^ 


discussed below, in which the effect of charge and Over¬ 
hauser noise on gate fidelities is (very nearly) additive. 

A simulation of the performance of our two-qubit gate 
in the presence of both charge and Overhauser noise is 
shown in figure In this figure we see that the trade¬ 
off of avoiding leakage at large while running the 

gate fast enough to avoid the accumulation of DC charge 
noise and Overhauser noise, leads to an optimal value for 
Jc/Jz of approximately 0.6 for large enough J^. Signi¬ 
ficantly, using just a single exchange pulse, we predict 
fidelities of ~ 0.95 with parameters currently accessible 
in experiment. 

This optimal value of Jd Jz is surprisingly large com¬ 
pared to intuitions garnered from prior arguments based 
on lowest order perturbation theory |28| . in which one 
must satisfy Jc ^ Jz- When this tighter contraint is sat¬ 
isfied, implementations of this operation would actually 
take longer than a more traditional approach involving 
multiple exchange pulses m-, and as seen in the results 
of our simulations, would in any case have low fidelities. 
That high fidelities are achievable with a large optimal 
value of Jc ~ 0.6 Jz is one of the main results of our study, 
and is due to the tuning of gate times and use of adiabatic 
pulses discussed in sections |IV A| and |IV B| respectively. 

There is still a large gap between the gate fidelit¬ 
ies found in the presence of these realistic noise sources 
and the upper limit imposed by leakage, as discussed in 
section |IV B| This could be mitigated in several ways. 
Firstly, one could increase our conservative choice of J^, 
with larger values of J^ leading to higher fidelities; with 
the caveat that it has been observed in experiment that 
qubit relaxation rates increase with Jz, likely due to 
phonons. More experimental investigation would be re¬ 
quired to determine the optimal choice of Jz, but ex- 


V. ALTERNATIVE GEOMETRIES 

While we have so far focussed mainly on the butterfly 
geometry due to its simplifying symmetries and higher 
performance, it seems likely that it will be much easier 
to fabricate a linear array of quantum dots. We there¬ 
fore extend our characterisation to include the anticip¬ 
ated performance of a linear geometry, as shown in figure 
121 

In section |IV B| we noted that the linear exchange 
coupling J 34 lacked the symmetry of the butterfly coup¬ 
ling, leading to an effective energy gap three times smal¬ 
ler than the butterfly geometry. Figure demonstrated 
that when this is taken into account in the absence of ex¬ 
perimental noise, and leakage is plotted as a function of 
Jc/AA, the leakage of the butterfly and linear arrange¬ 
ments is qualitatively similar. As detailed in section [llT| 
the Jj^ term in equation]^ is non-zero for the linear geo¬ 
metry, which makes it necessary to apply an echo pulse 
that anticommutes with the -\- cr^cr^ contribution 

j. Uz y y 

at t = 0 and at t = r/2 in order to perform a CPHASE 
gate (23]. One such choice (which we adopt in this work) 
is a pi rotation of the B qubit about the z-axis, which 
effects a operation. This could be experimentally 
implemented for the resonant exchange qubit by briefly 
shifting the value of Jzb- Since Jc changes the energy 
spectrum of the logical subspace, it is either necessary to 
ensure that Jc is turned off when the echo is applied. For 
our choice of a sinusoidal pulse envelope amounts to halv¬ 
ing the time of the pulse and repeating it twice. Prior to 
each of these Jc pulses one performs the single qubit azB 
operation. 

The smaller energy gap AA means that the gate must 
be run more slowly than in the butterfly geometry, mak¬ 
ing the linear geometry more susceptible to noise. The 
most significant contribution to infidelity in figure |5a] was 
DC noise on the intra-qubit couplings, which is geometry 
independent and diminishes fidelities for Jc/Jz < 0.5. 
Since leakage becomes significant in the linear geometry 
for Jc/Jz > 0.3, without suppressing DC noise, high- 
fidelity operation is not possible without modifying the 
gate. Here we demonstrate that suppression of DC noise 
is possible by modifying the echo pulse used. Applying a 
pulse at t = 0 and t = r /2 that anticommutes with the 
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Figure 8. (colour online) Entanglement fidelity for the lin¬ 
ear geometry subject to charge and Overhauser noise, as de¬ 
scribed in the main text, with noi se le vels calibrated to cor¬ 
respond to experiment (see table Till. A single a^ay echo 
pulse is applied at r/2 in order to echo out both intra-qubit 
psuedo-static charge noise and non two-qubit phase, 

which restores maximum fidelities to ~ 97%. Fidelities are 
limited predominantly by leakage (from above), Overhauser 
(from left and below), and high frequency noise (from below). 
The colourmap used is the same as in figure 


and af terms of equation]^ (such as a^a^) will echo 
out the effect of DC noise on the intra-qubit couplings. 
This can be merged with the af pulse discussed above. 
In this work we choose to apply a cr^a^ pulse, which be¬ 
comes (up to irrelevant phase) cr^(Jy when merged with 
af. Note that this amounts to an independent simul¬ 
taneous TT rotation for each qubit. In resonant exchange 
qubits this pulse can be achieved by oscillatory J pulses 
that are chosen to be out of phase by 7 r/ 2 . In the simula¬ 
tions shown in figure we demonstrate that these pulse 
sequences result in improved fidelities of ^ 97%. 

Similar to the butterfly geometry, gate performance is 
reduced for large and small values of Jc/Jz by leakage 
and (inter-qubit) DC charge noise respectively; and for 
small values of and/or Jd Jz by Overhauser noise. 

In this analysis, we have assumed that the echo pulses 
(which are single qubit operations) are instantaneous and 
performed with higher fidelity than two-qubit operations. 
In actual fact, current experiments demonstrate single 
qubit gate times several times longer in than the dur¬ 
ation of our two-qubit gate. At this preliminary stage, 
however, single qubit operations gate times have not been 
optimised, and we expect single qubit gate times could be 
reduced to a point where they are not limiting gate per¬ 
formance. In any case, our primary objective is to show 
that two-qubit gate operations in semiconductor systems 
can be implemented with times and fidelities comparable 
to single qubit operations, and so we have left a more 
detailed study to future work. 


VI. DISCUSSION 

In this paper we have extended the analysis of our pro¬ 
posal for a two-qubit CPHASE gate between resonant- 
exchange qubits m, which uses the large intra-qubit ex¬ 
change couplings to energetically suppress leakage caused 
by inter-qubit exchange couplings. In particular, we 
have: demonstrated that high-fidelity two-qubit opera¬ 
tions are possible with gate times comparable to single¬ 
qubit operations, provided that gate times are carefully 
tuned to account for perturbations beyond first order; 
demonstrated that leakage can be further suppressed by 
using adiabatic pulse profiles; and shown in conservat¬ 
ive simulations that the infidelities of our two-qubit gate 
are small (a few percent) even when subject to charge 
and Overhauser noise that has been calibrated to recent 
experiment. 

Our two-qubit gate works whenever there is an effective 
exchange coupling between the qubits. This could be a 
direct exchange coupling (as we have envisaged here), or 
an indirect coupling through an intermediate dot which 
has recently been shown to generate an effective exchange 
interaction 133110]. 

In our analysis, we found (intra-qubit) DC charge noise 
on the electrodes defining the quantum dots to be the 
most potent source of gate infidelity. While we have 
chosen our noise parameters to match experimental res¬ 
ults, we do not believe there is any fundamental reason 
why this noise cannot be substantially reduced. The 
next largest source of infidelity was the Overhauser field. 
We calibrated our simulations to reflect an unsuppressed 
Overhauser field in GaAs/AlGaAs heterostructures. It 
is possible to suppress the Overhauser field fluctuations 
using dynamic nuclear polarisation gUSSj, which can re¬ 
duce the standard deviation of the field fluctuations by 
several orders of magnitude (e.g. |3S|). It seems reas¬ 
onable to expect that our two-qubit gate’s fidelities may 
exceed 99% even without using echo pulses. Of course, 
since low-frequency (DC) noise and the Overhauser field 
fluctuations vary very slowly compared to gate times, it 
should be reasonably simple to use echo pulses or dynam¬ 
ical decoupling if necessary. 

A surprising result, perhaps, is that optimal ratio of 
inter-qubit coupling to intra-qubit coupling Jd Jz is ac¬ 
tually quite large, especially in the butterfly geometry. 
As a consequence, it is not necessary to run these gates 
particularly slowly, allowing our gate to run faster (and 
possibly with greater fidelity) than even the most care¬ 
fully constructed pulse sequences (e.g. [ 20 ]). 

Our approach of energetically suppressing spin-flip 
transitions in order to implement two-qubit gates using 
exchange coupling has utility in other qubit architectures, 
as we have already shown in singlet-triplet qubits |24j . 
The main benefit of this approach is to remove depend¬ 
ence on complicated pulse-sequences in order to achieve 
high fidelities. 

The relative simplicity of our two qubit gate, coupled 
with its high performance, commends it for implementa- 
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tion in contemporary experiment. 
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Supplementary A: Energy eigenstates for the rriz = 1 subspace 

In this section we list the energy eigenstates for the ruz = 1 subspace of the two resonant-exchange qubit system, in 
terms of the single qubit states listed in tableof the main text. Where an energy degeneracy occurs, we choose states 
that maximally exploit the symmetries of the butterfly geometry, in order to reproduce the splitting and coupling 
shown in figure of the main text. In particular, the states of the butterfly configuration are invariant under global 
spin rotations, allowing one to write the degenerate subspace in terms of states with definite angular momentum using 
standard Clebsch-Gordan coefficients. The states, along with their symmetries, are listed in the following table, in 
eigenstates are sorted by energy and then by total angular momentum S. Note that energies omit the constant —B± 
contribution shared by all states. The eigenstate label is chosen to indicate which elements of the single subspace 
are involved, and then append a colon followed by the total angular momentum (with the exception of logical states, 
where the total angular momentum is omitted). The “Parity” column indicates the parity accumulated by the state 
under the operation involving swapping dot 1 and 3 (the first sign) and 4 and 6 (the second sign) in the butterfly 
geometry. H—, for example, indicates that the state is unchanged under swapping the dots in the first qubit, but 
attracts a negative sign when swapping the dots of the second. The “Swap” column indicates the parity of the state 
when interchanging the qubits, where this is well-defined. Horizontal lines form the states into groups of equal energy. 


Label State Energy {+B±) S Parity Swap 
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Supplementary B: Derivation of the Effective Hamiltonian 

In this section we describe how the effective Hamiltonian of equations I and 2 in the main text is lifted from 
perturbation theory, the procedure for which follows standard practice |26| . The main idea is that, in the limit that 
inter-qubit coupling is small compared to intra-qubit coupling Jz, the evolution of the logical subspace in which 
we are interested should be well-approximated by low-order terms in a perturbation expansion around = 0; from 
which considerable insight might be gained into qubit dynamics. 

To generate the effective Hamiltonian for a given geometry, we perturb the Hamiltonian describing the two decoupled 
resonant exchange qubits with the Hamiltonian describing the two-qubit coupling of that geometry, to first order in 
a Rayleigh-Schroedinger perturbation expansion. We then construct the effective Hamiltonian using the resulting 
energies E and eigenvectors \E) that adiabatically map to the logical states using TLeff = This 

algorithm explicitly disregards any leakage operations, but captures the dominant dynamics on the logical subspace. 
The result is a diagonal Hamiltonian, except where the original logical states were degenerate. 

For our two resonant exchange qubit system, in which the |0I) and |10) states are degenerate (shown in figure 3 
of the main text), this results in a block-diagonal Hamiltonian when written in the basis {|00), |01) , |10) , |11)}. We 
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explicitly compute the effective Hamiltonians for the geometries shown in figure 2 of the main text; that is, the linear, 
butterfiy and rectangular geometries respectively. We also explicitly compute the coefficients of crfaf, , <Jy<7y, 

and a^af, which are then used to populate table II of the main text. 
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Note that to recover the provided coefficients in table II of the main text, we assume that all inter-qubit couplings 
(J 16 , J 2 b and J 34 ) are equal to Jc- 
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Figure 9. (colour online) In (a) we plot the ideal gate time 7!nax as computed using the “hybrid” method described in the text, 
which is our best approach for maximising the fidelity of gate operations. 71nax(Jc/Jz) is defined as in the text, and plotted 
here for = 1.65/rel/. In (b) we plot the relative difference T/7(nax of the other methods described in the text, namely 
“adiabatic” (green dashed) and “first order” (red dotted) methods, to the “hybrid” estimate (blue solid). The plots share the 
same x axis, which is the ratio of intra-qubit coupling to inter-qubit coupling Jc/Jz. Note that the first order approximation 
consistently underestimates the phase accrual rate, and so overestimates the gate time. While an adiabatic phase estimation 
does better, it ends up underestimating the phase due to the importance of evolution off the logical subspace for large enough 

Jc 


Supplementary C: Estimating Ideal Gate Time 

As mentioned in section IVA, over- or under-estimating the ideal gate time r leads to over- or under-accrual of 
two-qubit phase, resulting in poor gate fidelities. In particular, using a first order approximation neglects higher order 
terms that give rise to convex non-linearity in Jzz(Jc), and so consistently exaggerates gate times. In this section, we 
describe the method by which we obtain accurate timing estimates for use in the simulations of the main text. 

There are two methods that one might employ to find an estimate of t that includes high order perturbations. 
One can make the assumption that evolution is perfectly adiabatic, solve for the eigenvalues at different values of 
inter-qubit coupling Jc, and thus infer the two-qubit phase accumulation rate Jzzl or, alternatively, one can locally 
maximise the gate fidelity over gate times. The first method works well for small values of Jd Jz (where the adiabatic 
approximation make sense), and the second method works well for larger values of Jd Jz (where gate times are shorter, 
and numerical integration has less time to accumulate error). We therefore use a hybrid approach that works well 
across all values of Jd Jz'- for Jd Jz < 0.15, we use the adiabatic approach, and for Jd Jz > 0.15 we use the adiabatic 
approach to seed a numerical optimisation of gate fidelity over gate times. 

As noted in the main text, the non-linearity of Jzz{Jc) means that we must repeat this estimation process for each 
pulse shape considered. As the dynamics of our (noiseless) two-qubit system is determined by the ratio of Jc/Jz, 
given a particular pulse shape, we need only optimise over the ratio of the two physical degrees of freedom, rather 
than both. We therefore proxy the optimisation of t( Jc, Jz) for any given pulse shape by optimisations over the single 
parameter function T{JdJz)- In practice we fix a value of Jz = J™^ and vary J^- t can be recovered from T using: 

Tref 

T{Jc,Jz) = ^r{JdJz)- 

^ z 

The gate time for maximum fidelity 7)nax for the butterfly configuration using a sinusoidal adiabatic pulse profile is 
shown in figure along with the relative error of the first order and adiabatic estimations. We note that the first 
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order approximation consistently underestimates the phase accrual rate, and so overestimates the gate time. While 
an adiabatic phase estimation does better at first, it ends up underestimating the phase due to the importance of 
evolution off the logical subspace for large enough J^- 


Supplementary D: Accelerating Monte-Carlo Convergence in Overhauser Simulations 


Simulations involving the Overhauser field modelled as described in the main text requires averaging over a multi¬ 
variate normal distribution in 6 variables (11 if DC charge noise is also considered). While a naive Monte-Carlo 
simulation that samples the local Overhauser field contribution for each dot will converge eventually, we speed up 
the convergence by instead sampling preferentially from linear combinations of local Overhauser contributions that 
appear at lower order in perturbation theory. 

Inspired by perturbation theory, we form a spanning basis for the magnetic fields: 



n 


yn<3 n>4 J 
Aij = Bi — Bj for A 13 and A 46 
Aij]^ — Bj^ ‘^Bj -t- Bf^ for A 123 and A 455 


where Bn (with n € [1,6]) is the z-component of the Overhauser field at dot n. We learn from perturbation theory 
that A 123 and A 456 first contribute at zeroth order in perturbation theory, A 13 and A 46 at first order, and AB at 
fourth order. As Bq is a global field, it does not contribute. 

The higher the order at which the terms contribute, the more significant that term is to the dynamics of the qubit 
system. It makes sense, therefore, to sample lower order terms more often. Technically, one should weight each term 
roughly as {Jc/Jz)°, where o is the order at which the term appears in perturbation theory. We found it simpler, 
however, to use a conservative ratio of 5 between terms of different order. That is, for each 100 simulations, we sample 
Ai 23 and A 456 100 times, A 13 and A 46 20 times, and AB 4 times. 

As the terms of the new basis are linear compositions of the Bn terms, the random variable associated with each 
term will differ from the underlying local field fluctuations. In particular, if each Bn ^ JV{Bq, a^), then the relevant 
terms of new basis will be sampled from the distributions: 

Aijk ^Af{0,4a%) 

A,, ~Af( 0 , 24 ) 

AB ~Af(0,2(T|) 


Supplementary E: Cross-Coupling 


Perturbation theory provides the insight that beyond first order the Overhauser field cross-couples with charge noise 
in Jc at second order and above. This allows the Overhauser field to contribute non-trivially to two-qubit phases. For 
example, the two-qubit phase accumulation rate to second order for the butterfly configuration is given by: 


+ 243J/Jf (J/ + ) 

j f ^123 , A 456 

* 81''^ ^ 

+0( JeA^) + 0{J^A) + 0{Jl), 


where Aijk = Bi — 2Bj + Bk, and where the coloured and starred line is the effect of cross-coupling. The intuition is 
that as Jc/ Jz and the as increase, the greater the two-qubit cross-coupling error. For the useful range of Jd Jz values, 
and the experimentally inspired values oi as used in the main text, the effect of this cross-coupling is negligible. 





